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Abstract:  In  the  work,  transmission  of  nerve  impulses  along  nerve  fibers  is  simulated.  Research  is  being  conducted  in 
the  framework  of  the  electrical  theory  of  propagation  of  the  action  potential.  The  soliton  approach  is  used.  The  numerical 
experiment  on  the  evolution  of  the  nerve  impulse  and  the  laws  of  motion  is  conducted.  It  is  shown  that  the  stable  form  of 
nerve  impulses  solitons  is  realized  under  different  initial  conditions.  The  threshold  character  of  occurrence  of  nerve  pulse  is 
simulated.  It  is  shown  that  the  number  of  pulses  produced  changes  depending  on  the  degree  of  nonlinearity:  with  strong 
initial  excitation  of  large  amplitude  soliton  is  unstable;  it  breaks  up  into  a multitude  number  of  solitons  with  small  ampli- 
tudes. In  this  case,  the  greater  the  nonlinearity  parameter,  the  greater  the  number  of  births  of  solitons.  Unusual  solitonlike 
regimes  of  interaction  of  nonlinear  pulses  excitations  are  illustrated;  in  certain  anomalies,  colliding  nerve  impulses  are  re- 
flected instead  of  their  usual  quenching.  The  possibility  of  the  decay  of  the  nerve  impulse  at  the  bifurcation  of  the  nerve 
fibers  or  the  appearance  of  heterogeneity  of  the  passing  of  a nerve  impulse  (the  presence  of  dissipation  in  the  environment) 
is  investigated.  Physical  modeling  allows  obtaining  and  studying  the  analytical  results  and  elucidating  the  physical  prin- 
ciples of  biological  processes. 
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1.  Introduction 

To  coordinate  work  of  the  different  systems  in  the  hu- 
man body  requires  harmonious  work  of  its  most  complex 
system  - the  nervous  system.  Its  basic  function  is  the  prop- 
agation of  an  electric  pulse  through  the  nerve  fiber.  Physio- 
logical studies  have  shown  that  propagation  of  a nervous 
pulse  through  neural  nets  has  auto  wave  character.  T. 
Heimburg  and  A.D.  Jackson  proposed  a new  neural  theory 
called  Soliton  theory.  In  this  theory,  the  nerve  conduction 
is  proposed  as  a density  wave.  [1]  B.Lautrup,  R.Appali 
and  others  [2]  examined  a class  of  solitons,  obtained  from 
a generalization  of  the  Boussinesq  equation.  Soliton  model 
is  based  on  the  propagation  of  a localized  density  wave  in 
the  axon  membrane.  This  allows  simulating  such  processes 
using  Soliton  Theory,  which  describes  the  form  and  veloci- 
ty of  propagation  of  the  nervous  pulse,  stability  of  its  cha- 
racteristic parameters  in  time,  and  process  locality  with 
perfect  precision.  [3] 

According  to  electrical  theory,  the  process  of  transmis- 
sion of  the  pulse  along  the  nerve  fiber  is  a process  of  prop- 
agating the  action  potential:  the  nerve  fiber  extends  the 
electrochemical  reaction,  which  generates  a running  im- 
pulse of  voltage.  The  nerve  pulse  or  action  potential  propa- 


gates through  the  nerve  fiber  without  decay  and  with  con- 
stant velocity.  An  impulse  can  be  formed  and  spread  be- 
cause the  nerve  fiber  contains  a nonlinear  element,  which 
suppresses  small  deviations  from  normal  state  and  streng- 
thens big  ones.  In  this  case,  the  nonlinear  dependence  of 
the  membrane  permeability  (through  which  the  nerve  im- 
pulse is  spread)  from  the  momentum  and  the  diffusion  of 
ions  across  the  membrane  are  balanced.  In  the  paper  [4]  is 
shown  that  nonlinearity  can  actually  balance  diffusion  and, 
as  a result,  a running  solitary  wave  with  constant  speed  and 
shape  can  occur,  i.e.  a soliton.  Thus  a nerve  impulse  is  an 
extension  of  an  action  potential.  Physiological  studies  have 
shown  that  propagation  of  a nervous  pulse  through  neural 
nets  has  autowave  character.  This  also  allows  simulating 
such  processes  using  Soliton  Theory  which  describes  the 
form  and  velocity  of  propagation  of  the  nervous  pulse,  sta- 
bility of  its  characteristic  parameters  in  time,  and  process 
locality  with  very  good  precision.  [5] 

2.  Theory  and  Numerical  Experiment 

Since  the  nerve  impulse  is  resistant,  it  extends  without 
distortions  and  significant  loss  of  energy;  its  movements 
can  be  modeled  as  a Korteweg  de  Vries  soliton.  In  this  case 
it  is  described  by  the  nonlinear  differential  equation  [6]: 
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du  d3u  du 


0) 


u(x,  0)  = u0  exp 


(3) 


where  the  second  term  characterizes  a weak  (cubic)  dis- 
persion, and  the  third  term  describes  nonlinearity.  The 
solution  of  this  equation  has  the  following  form, 


U2(X~Vt\ 

ch 

where  u is  the  action  potential , umax  is  the  amplitude  of 
the  action  potential  (in  our  case  of  the  nerve  pulse),  v is  the 
velocity 


(v  Umax/ 3)  (2  Cl) 

, and  A is  determined  from  the  ratio 

A=  V 12/3  Aw,  (2  b). 

From  (2  a)  it  follows,  that  the  propagation  velocity  of  a 
pulse  increases  with  increase  in  the  amplitude,  whereas 
from  (2b)  it  follows  that  the  soliton  width  decreases  with 
the  increase  in  the  amplitude.  Therefore,  the  solitons  line 
up  in  the  order  of  decreasing  amplitudes  and  growing 
widths  “Fig.  1”. 


Figure  1.  Soliton  distributions  at  different  times.  (Numerical  experiment). 

We  performed  a numerical  experiment  of  formation, 
evolution  and  movement  mechanism  of  solitons.  The  stable 
form  of  solitons  can  be  obtained  with  different  initial  con- 
ditions. 

We  chose  the  initial  perturbation  in  two  forms: 


(where  / is  the  linear  size  of  initial  perturbation)  and  a 
hyperbolic  cosine.  However,  regardless  of  the  character  of 
initial  perturbation,  over  time  the  pulses  transform  into  a 
soliton.  Obviously,  the  initial  conditions  are  connected 
with  the  signal  received  from  the  brain  through  conducting 
ways  and  depend  on  the  state  of  the  organs  or  the  mental 
state  of  the  person.  ( These  data  are  consistent  with  [6] ) 
English  physiologists  A.  Hodgkin  and  A.  Huxley  have 
developed  a theory  (based  on  numerous  experiments) 
where  they  concluded  that  the  shape  and  the  speed  of  an 
impulse  don’t  depend  on  the  stimulus  of  the  nerve.  If  the 
stimulus  is  high,  several  consecutive  impulses  are  produced, 
while  if  it  is  low,  no  impulse  will  be  produced  through  the 
nerve;  the  minimal  power  of  stimulus  to  produce  an 
impulse  is  called  threshold.  Each  impulse  transfers  one  unit 
of  information  and  the  receiving  devices  count  the  number 
of  elementary  particles  of  information  received  and  the 
time  it  took  [7].  All  this  reminds  of  the  movement  of 
solitons.  In  paper  [8]  Hodgkin-Huxley  and  Soliton  theories 
are  described  and  a theoretical  comparison  has  been  carried 
out  throughout  the  analysis  of  the  theories  and  models. 

In  the  Korteweg  de  Vries  equation  the  degree  of  nonli- 
nearity is  characterized  by  a dimensionless  parameter 
a = l(u0/ fl)1/ 2 . The  numerical  experiment  showed  that  at 
a < (Jmiri  the  solitons  are  absent, and  instead  sine  wave 
oscillations  arise  “Fig.  2”. 


UL 


Figure  2.  Sine  wave  oscillations. o < omil 
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From  “Fig  2a”  it  is  visible  that  solitons  are  absent,  but  to 
the  left  of  an  initial  impulse  there  are  oscillations.  This  case 
is  realized  at  /3  > 0.  At  /3  < 0 (“Fig  2b”)  oscillations 
arise  on  the  right.  The  received  result  follows  from  Korte- 
weg  de  Vries  equation. 

We  investigate  character  of  the  solution  of  the  equation 
(1)  depending  on  values  a.  For  a = 5.9  we  studied  evolu- 
tion of  an  initial  impulse  for  times  tx  < t2  < t3 . Thus  the 
initial  impulse  moves  to  the  right,  the  first  solitons  ampli- 
tude grows  over  time,  at  t2  the  second  soliton  is  born 
whose  amplitude  also  grows  over  time;  at  t3  we  have  2 
solitons,  and  the  amplitude  of  the  first  one  is  higher  than 
the  amplitude  of  the  second  one.  At  a = 9 and  the  same 
initial  condition  the  decision  for  the  same  values  of  time  is 
received.  It  can  be  seen  from  the  figure  that  2 separate  soli- 
tons on  the  right  and  the  third,  a new  one,  on  the  left  are 
visible.  Thus,  with  increase  of  a at  an  identical  initial  im- 
pulse the  number  of  born  solitons  increases  “Fig.  3”. 

It  is  known  that  in  vivo  (living  organism)  nerve  impulses 
moving  along  the  axon  toward  each  other  collide  and  anni- 
hilate each  other,  which  provides  one-way  transmission  of 
signals  from  receptors  in  the  central  nervous  system.  How- 
ever, the  works  of  Aslamanidi  O.V.,  Momev  O.A.  [9,10] 
are  based  on  numerical  experiments  with  Hodgkin-Huxley 
equations  predict  a new  effect:  for  some  values  of  the  pa- 
rameters colliding  nerve  impulses  may  reflect.  It  is  shown 
that  in  the  soliton  mode,  the  propagating  impulse  is  a 
“doublet”  consisting  of  a high-amplitude  impulse-leader 
and  an  accompanying  low-amplitude  wave.  The  first  colli- 
sion of  the  leaders  causes  annihilation,  and  the  second  col- 
lision of  the  low-amplitude  waves  leads  to  summation.  As 
a result,  at  the  point  of  collision  the  potential  reaches  the 
above-threshold  value,  causing  regeneration  of  “doublets” 
that  become  scattered  in  different  directions. 
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Figure  4.  Collision  of  two  solitons.  The  figures  show  successive  moments 
of  time. 


Figure  3.  The  different  number  of  originating  solitons  for  different  values 
of  nonlinear  parameter. 

We  have  modeled  the  collision  of  two  solitons  occurring 
in  a similar  mode,  on  the  basis  of  Korteweg  de  Vries  equa- 
tion “Fig.  4”. 


u 


The  analysis  of  such  regimes  of  interaction  of  nonlinear 
pulse  can  be  useful  in  deciphering  the  mechanisms  of  car- 
diac arrhythmias.  It  is  known  that  the  Korteweg  de  Vries 
soliton  gradually  slows  down  under  friction,  reducing 
height  and  becoming  blurred.  The  lifetime  is  inversely  pro- 
portional to  the  force  of  friction.  Analogically,  when  the 
cross  section  is  changed  or  nerve  fibers  are  bifurcated,  the 
passage  of  nerve  impulses  may  become  difficult  or  even 
completely  blocked.  In  the  equation,  this  can  be  considered 
as  the  presence  of  friction.  [6]  We  have  also  investigated 
the  solitons’  evolution  in  the  presence  of  friction.  When 
taking  friction  into  account,  the  Korteweg  de  Vries  equa- 
tion becomes  the  form 


du  du  d2u  d3u 
at  +U5t  ~U'dtI  + ^dx^  = ° 


ti 


u 


t2 


u 


where  ft  is  the  coefficient  of  friction.  Suppose  that  at 
t = 0 we  have  an  initial  perturbation  in  the  form 

u(x,  0)  = u0  ch_2(x/0- 

At  (i  = 2 the  numerical  experiment  was  made  for  t=  tl< 
t2<  t3<  t4<  t5.  “Fig. 5a”.  Next,  the  same  experiment  was 
performed  for  (i  = 10  and  for  the  same  values  of  time.  As 
expected,  attenuation  occurs  faster  with  with  the  increasing 
coefficient  of  friction.  “Fig.5b”. 
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Figure  5.  Solitons  evolution  at  different  coefficients  of  friction  ju  for  the 
same  values  of  time. 

3.  Conclusion 

Thus,  different  regimes  of  propagation  of  nerve  pulses 
can  be  realized  under  different  conditions.  Solitonlike  solu- 
tions describe  a wide  range  of  phenomena  that  accompany 
the  transmission  of  nerve  pulses,  heart  rate,  self-oscillation. 
The  soliton  models  of  processes  make  it  possible  to  obtain 
adequate  description  and  identification  of  the  state  of  such 
systems  as  the  conducting  muscular  and  circulatory  sys- 
tems of  the  heart.  [11]  By  changing  parameters  of  the  soli- 
ton model,  it  is  possible  to  model  the  pathologies  of  nerv- 
ous conducting  ways  of  the  vascular  and  muscular  system 
of  the  heart.  The  soliton  approach  opens  up  large  possibili- 
ties in  the  research  into  the  complex  processes  in  living 
organisms.  Possible  pathology  in  the  vascular  system 
changes  spectral  picture  of  the  soliton  waves  and  allows 
diagnosing  a variety  of  diseases  of  nerve  pathways  and  the 
circulatory  system.  Mathematical  modeling  provides  the 
analytical  results  and  the  way  to  identify  the  physical  prin- 
ciples of  biological  processes. 
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